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3 ABSTRACT 

4 Studies over a period of several decades have resulted in a relatively simple set of equations 

5 describing the tidally and width-averaged balances of momentum and salt in a rectangular 

6 estuary. We rewrite these equations in a fully non-dimensional form that yields two non- 

7 dimensional variables: (i) the estuarine Froude number; and (ii) a modified tidal Froude 

8 number. The latter is the product of the tidal Froude number and the square root of 

9 the estuarine aspect ratio. These two variables are used to define a prognostic estuary 
10 classification scheme, which compares favourably with published estuarine data. 



1 1. Introduction 



2 Since the introduction of stratification-circulation diagram by Hansen and Rattray ( 1966 ) , 

3 numerous estuarine classification schemes have been proposed. The reader might ask - why 

4 revisit this topic? Our motivation for pursuing a new classification scheme stems from no- 



table recent advances in estuarine physics, many of which are reviewed in MacCready and 



6 [Geyer (2010). These advances led us to hypothesize that there might be a simple means to 



7 determine the conditions under which a sufficiently well behaved estuary will be well mixed, 

8 partially mixed, or highly stratified. We start by outlining the classical tidally averaged 



9 model as presented by MacCready and Geyer (2010). We then rewrite the equations of this 



20 model in non-dimensional form. Using this new set of equations we develop our classification 

21 scheme, and then compare its predictions with field observations. 

22 2. Classical Tidally Averaged Model 

23 The physics of estuarine circulation is governed by the competing infiuences of river 

24 and oceanic fiows. While the former adds fresh water, the latter adds denser salt water 

25 which moves landward due to the combined effect of tides and gravitational circulation (or 

26 exchange flow). The complicated balance between the river, the exchange flow and the tides 

27 determines the estuarine velocity and salinity structure. 

28 We consider an idealized rectangular estuary of depth H and width B. The origin of 

29 the coordinate system is at the free surface at the mouth of the estuary with the horizontal 

30 {x) axis pointing seawards and the vertical {z) axis pointing upwards. Therefore, both the 

31 horizontal and vertical distances within the estuary are negative quantities. To obtain the 

32 width-averaged and tidally-averaged horizontal velocity (u), and salinity (s) distribution in 

33 the estuary, these quantities are flrst decomposed into depth averaged (overbar) and depth 

34 varying (prime) components: u = u{x,t) + u'{x,z,t), s = s{x,t) + s'{x,z,t). The quantity 

35 M = Qr/A is the cross-sectionally averaged river velocity, where Qr is the mean river flow 



36 rate and A = BH . The solution for both partial and well mixed estuaries was given by 

37 Hansen and Rattray (1965) (for recent review, see MacCready and Geyer ( 2010[ )): 



U = U + U = uPi + UeP2 



(1) 

(2) 
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In (pi), ^ = z/H G [—1,0] is the normalized vertical coordinate. The subscript x implies 
d/dx where x is dimensional. Ks is the vertical eddy diffusivity. For exchange dominated 
estuaries, an important parameter is the exchange velocity scale: 



ue = c'H'J:J {A8K_ 



M 



(4) 



49 Here c = yJg[3socnPl^ is twice the speed of the fastest internal wave that can be supported in 

50 an estuary ( MacCready and Geyer|2010 ). Km is the vertical eddy viscosity and /3 = 7.7 x 10"^ 

51 psu^^. The non-dimensional salinity is defined as S = s/socn, where Socn is the ocean salinity. 

52 Equations ((l])-([2]) were derived under the assumption that the density field is governed 

53 by the linear equation of state: p = po{l + (3s) where po is the density of fresh water. The 



54 details of the derivation are well documented in MacCready (1999, 2004) 

55 The salt balance is given by: 

d 



dt 



Tidx = — u'T,' + Kh^x ~uTi 

■' 

accumulation exchange tidal river 



(5) 



58 where, Kh is the horizontal diffusivity. This equation physically implies that the temporal 

59 salt accumulation in an estuary is due to the competition between salt addition and removal 



60 processes. While exchange (note that u'YJ is negative) and tidal processes add salt, river 

61 inflow removes it. At steady state (IS]) can be rewritten as: 



63 
64 
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where 



Lh = Kh/u 

Lei = 0.0l9uH^/Ks 

LE2 = 0.031cHy{KsKMy^' 



Le3 = 0.024{c/uY/'cHy [KsKl] 
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(7) 



69 The different terms in ([6]) are as follows: R is the river term, T is the tidal term, while Ei^ 



E2 and E3 are the different components of the exchange term. Hansen and Rattray (1965) 



71 presented ([6]) in a slightly different form, and MacCready (2004, 2007) introduced the length 

72 scales in ([T]). 

73 The length scales in ([T]) depend upon the mixing co-efficients: Ks, Km and Kh- Making 



74 use of an extensive study of Willapa Bay, Banas et al. (2004) proposed: 



Kh = aiUrB; 



76 where Oi = 0.035 and ut is the amplitude of the depth averaged tidal flow. Based on field 



77 studies and modeling of the Hudson River estuary, Ralston et al. (2008) obtained 



Km = aoCoUTH and Ks = Km/Sc; 



(9) 



79 where ao = 0.028, Cd = 0.0026 and Sc = 2.2 is a Schmidt number. We will use (g and (g 

80 in the development of a non-dimensional set of equations. 

81 While the governing equations ([T|, (|2| and ^ are elegant representations of the prob- 

82 lem of estuarine circulation, they are sufficiently complicated that simplifications have been 



83 sought after. Numerous investigators, including Hansen and Rattray (1965); ?); ?); Mac- 



Cready (2004); MacCready and Geyer (2010) have assumed u <C ue, which yields: 



S 
R 



[Le3^x} + Lh^3 
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(10) 



87 ? further reduced (10) to two simple cases with analytical solutions, the exchange dominated 



88 case (T — 7- 0), and the tidally dominated case {E^ — )■ 0). While these approximations have 

89 been widely used there does not appear to have been any serious attempt to determine the 

90 conditions under which they are applicable. 



91 3. Non-dimensional Tidally Averaged Model 

92 In this section we rewrite the governing equations M, d2]) and (^ in non-dimensional 

93 form in anticipation of (i) revealing the important non-dimensional parameters governing 

94 the problem, and (ii) facilitating comparison of the relative magnitude of each of the terms 

95 in ([6]). Defining X = x/Les, (|6| can be rewritten as: 



96 
97 
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(12) 



Le3 V 0-024 
101 The velocity u and ut have been non-dimensionalized by c to obtain the densimetric estuarine 



Froude number Fr = u/c and the tidal Froude number Ft = ut/c. Substituting (12) into 



(11) yields: 
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106 where Ci = 2.17, C2 = 1.34, C3 = 8.16 x 10~^, and the modified tidal Froude number, 

107 Ft = Ft\/B/H. Typically the estuarine aspect ratio B/H ~ 0(10^ — 10^), see Table 

108 fl| The magnitude of different terms in ( p!3| can be easily compared by noting that < 

109 O {Fr^/^) < O [Fr"^/^) < O (1) < O (Fr^^/^). The tidal term (T) however depends on an 
no additional parameter Ft-, whose (order of) magnitude needs to be known for making the 

111 comparison. 

112 Like the salt balance equation, the momentum and salinity equations, i.e. ([I]) and ([2]) 

113 can also be expressed in non-dimensional form as follows: 



115 
116 



U — CiFp^ T,x P2 

S = S + C^F^ Ej^ P4 + C^Fp^ Sx P3 
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(14) 
(15) 



The constants C4 = 0.667, C5 = 47.0 and Cg = 70.5. In (14), the quantity U = u/c is the 
non-dimensional horizontal velocity (not to be confused with Fr, which is u/c). Equations 



119 (13)-(15) are the non-dimensional governing equations for our idealized estuary. 



Eq. (13) poses a non-linear initial value problem which can only be solved numerically. 



For that, the conditions at the estuary mouth have to be determined. One such condition 
is S (0, —1) = 1; meaning the salinity at the bed of the estuary at its mouth has to be the 



same as the ocean salinity. Substituting (15) into (13) and making use of this condition, we 
obtain 

{ExloY + CrF'J' (Sxlo)' + (CsF^/' + CsFr' F-'^') S^lo = 1; (16) 



126 where C7 = 5.31 and Cg = 6.04. Eq. (16) is actually the non-dimensional version of Eq. (19) 



127 of MacCready (2004). Being a cubic equation, it can be solved analytically to evaluate the 



128 salinity gradient at the estuary mouth, Sx|o- Additionally, (16) indicates that Sx|o is only 



a function of Fr and Ft- The variation of Sx|o with these two Froude numbers is depicted 
in Fig. QJ The figure shows that < T.x\o < 1 over the entire parameter space. 



4. Estuary Classification 



132 Our goal is to develop a simple classification scheme that distinguishes between well- 

133 mixed, partially mixed and highly stratified estuaries. A relevant parameter for classifying 

134 estuaries is the non-dimensional salinity stratification at the estuary mouth, $o- It is defined 

135 as follows: 

136 $o = S(0,-l)-S(0,0) (17) 

137 This parameter ranges between and 1. While the lower limit implies a very well mixed 



estuary, the upper limit indicates the transition to salt wedge. Substituting (15) into (17) 
yields: 

$0 = CgF^/' (Sxlo)' + CioF^/'SxIo; (18) 

where Cg = 7.06 and Cio = 8.82. If Fr and Ft are known, then Sx|o can be directly 



142 obtained by solving (16). Consequently, $o can be evaluated from (18), yielding Fig. [2j 

143 We follow Hansen and Rattray (1966) and use the condition $o = 0.1 to define the 



transition between well mixed and partially mixed estuaries. To distinguish between partially 
mixed and highly stratified estuaries we use the condition $0 = 1-0, corresponding to fresh 
surface water extending to the mouth of the estuary. Our classification scheme is obtained 
by plotting these transitional criteria on Fig. |2| When F^ = 0, the transition between 
well-mixed and partially-mixed estuaries is predicted to occur at Fji = 0.0017, and from 
partially-mixed to highly stratified at F^ = 0.113. The value of Fr for both transitions 
increases as Ft increases, the increase being more rapid for the transition from partially- 
mixed to highly stratified estuaries. These results are in qualitative agreement with Fig. 2.7 



152 of Geyer (2010). 



5. Discussion 



154 Together (16) and (18) provide new insight into estuarine physics. Apart from broadly 



classifying estuaries into three categories, viz. highly stratified, partially mixed and well 
mixed, the equation set identifies Fr and Ft to be the only two parameters determining 
the stratification at the estuary mouth, $o- The new non-dimensional parameter Ft = 
Ft\/B /H reveals that "tidal effect" is not simply represented by the tidal Froude number Ft-, 
but the latter combined with the square-root of the estuarine aspect ratio B/H. Moreover 
the equation set predicts $o; given Fr and Ft- If estuarine condition changes, e.g. river flow 
changes from low to high, tidal flow changes from spring to neap, or estuary depth changes 
due to dredging, the parameters Fr and Ft will change correspondingly. These newly 
obtained Froude numbers will produce a new $o, which reflects the response of estuarine 
circulation and mixing to variability. 

To test the applicability of our classification scheme we made use of the field data pre- 
sented in 



Prandle 



(1985). Using these data we have computed Fr, Ft, B/H, and Ft directly, 

167 and $0 from (16) and (18); see Table [I] We have compared the computed value of $0 with 

168 the measured value in Fig. [3j The comparison is good considering the accuracy to which Fr 

169 and Ft can be determined from field data. 



It is interesting to note that if both Fr and Ft are small then (16) reduces to Sx|o = 1 



171 and (18) reduces to: 



$r 



7F, 



2/3. 



R ' 



(19) 



which is exactly the same as Eq. (19) of MacCready and Geyer (2010). They have found 
this equation by combining Knudsen's relations (Knudsen 1900) with (J9]). Since Knudsen's 
relations are derived from mass and salt balances and do not consider momentum balance, 



176 (19) provides a rather simplistic prediction. We compare (19) with the exact solution of 



(18) for both Ft = and Ft = 30 in Fig. 4 When Ft = 0, (19) is very accurate up to 



178 the transition between well- mixed and partially mixed estuaries ($0 = 0.1). The deviation 



179 between (19) and the exact solution increases with increasing Fr. Eq. (19) predicts the 



transition between partially mixed and highly stratified estuaries ($o = 1-0) at Fr = 0.054 



rather than at Fr = 0.113. When Ft = 30, (19) is always less than the exact solution in 
terms of $o- 



183 We also compare the theoretical results with the field data of Prandle (1985) in Table 



[I] and Fig. 2.6 of Geyer (2010) in Fig. Ill This comparison is mainly intended to provide a 



185 qualitative estimate. We have chosen to plot (18) for Ft = and 30, since estuaries mostly 



have Ft within this range. Ideally, most of the partially and well mixed estuaries should 
cluster within the grey region bounded by the lines Ft = and 30, which is indeed the 
case. The most important aspect of this comparison is that the theoretical curves follow 
the overall trend of the field data. However these curves grossly over-predict $o, therefore 
they under-predict vertical mixing. This discrepancy may arise if the values of $o were 
measured at an upstream location, rather than at the mouth (which might be the case for 



192 the data points of Geyer (2010)). Surprisingly, disagreement between theory and field data 



193 does not appear in Fig. 2.6 of Geyer (2010). The latter figure shows that the line, referred 



194 to as "Eq. (2.22)", matches very well with the data points. Although Eq. (2.22) is actually 



$0 = ^■'^'^Pr (we calculated the coefficient from the associated text in 



Geyer (2010)), it is 



mistakenly plotted as $o ~ 3Fr 



R 



Finally we refer to the assumptions behind our theoretical analyses and their conse- 
quences. We have simplified the problem by assuming a tidally averaged estuary with rect- 
angular geometry. In real estuaries bathymetry can play a crucial role in determining the 
estuarine circulation. Moreover the appearance of just two parameters {Fr and Ft) in our 
equations is a consequence of the empirical equations ([s]) and ^. These two equations 
are also used in determining the coefficients Ci, C2, . . . , Cio- All these coefficients are found 
to depend upon Sc, making it the most important parameter in this regard; see Table [2] 
Following Ralston et al. (2008), Sc = 2.2 in all our calculations. Although ([s]) and ^ 
are simple and elegant, they may not be very realistic. In real estuaries both Km and Ks 



206 are variables. Moreover, other empirical parameterizations have shown that Km depends 

207 upon Richardson number (MacCready and GeyerpOlO). While the inclusion of any relevant 



208 third parameter might improve the predictability of the classification scheme, the value of 

209 this improvement would have to be weighed against the added complexity of the resulting 

210 classification scheme. 

211 6. Conclusions 

212 The equations governing the physics of estuarine circulation have been presented in non- 
213 dimensional form. The two resulting non-dimensional parameters are the estuarine Froude 
214 number, Fr, and the modified tidal Froude number. Ft- Given these parameters the non- 
215 dimensional salinity gradient at the estuary mouth, Sxlo; and the non-dimensional salinity 

216 stratification (also at the estuary mouth), $0, can be computed. The latter result forms the 

217 basis of a classification scheme that can be used to predict whether an estuary is fully or 

218 partially mixed, or highly stratified. The predictions of this classification scheme compare 

219 well with estuarine data. 
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Table 1. Estimates of estuarine parameters calculated using the data of Prandle (1985) 



and values of B obtained from maps. Eq. (18) is used to obtain $o-theory. 



Estuary Name 



Ft B/H Ft $o $o-theory 



Vellar 


1.27 


0.64 


200 


9.0 


1.00 


1.00 


Columbia 


0.026 


0.43 


150 


5.3 


0.40 


0.50 


James 


0.004 


0.25 


360 


4.7 


0.22 


0.17 


Tees 


0.014 


1.03 


75 


8.9 


0.18 


0.33 


Southampton Waterway 


0.0012 


0.37 


200 


5.2 


0.10 


0.06 


Tay 


0.014 


1.38 


400 


27 


0.10 


0.17 


Narrows of the Mersey 


0.0009 


0.83 


65 


6.7 


0.05 


0.05 


Bristol Channel 


0.006 


1.59 


300 


27 


0.02 


0.06 
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Table 2. List of coefficients used in different equations. 



Coefficient 


Value 


c. 


imsc^i^ 


C2 


0.7925c2/3 


Cs 


41.7aoaiCfl5c-i/3 


c. 


0.86850-1/=^ 


c. 


36.2Sci/3 


c. 


41.7^c2/3 


Cj 


4.085ci/3 


Cs 


3.575c2/3 


c. 


5.43^c^/3 


Cio 


5.215c2/3 
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List of Figures 



The variation of salinity gradient at the estuary mouth (Sx |o) with the estuary 

control variables - Fr and Ft- The solid lines represent isocontours of Sx|o- 

Estuary Classification Diagram. The lines represent isocontours of $o- The 

three regions represent three types of estuaries: (a) light Grey- well mixed, (b) 

White - partially mixed, and (c) dark Grey- highly stratified or salt wedge. 

The letters denote estuaries: C - Columbia, J - James, Te - Tees, SW - 

Southampton Waterway, Ta - Tay, NM - Narrows of the Mersey and B - 

Bristol Channel. For data, see Table [1} 

Comparison between stratification at the estuary mouth obtained from theory 

with field data). 

Comparison between our estuary classification scheme and the approximation 



.2/3 



$0 = 7Ffj in (19). The grey area indicates the region where estuaries should 



ideally cluster. Field data from Geyer (2010) and Prandle (1985) are plotted 



for comparison with the theoretical predictions. 



HZ] 



m 



m 



16 



30 



t:^20 



^10 







y ■ ■ ■ 1 
/O.I 


/ 


/ / 


y\ • 1 r 




/ 


/ 

0.3 

/ 


0.5 
/ 


^0.7 \ 




0.3 




^ 




1 0.5 




.... 1 


^ 


0.9--. 


"^ ' 


. 1 







10 



-3 



10 



-2 



F 



R 



10 



-1 



Fig. 1. The variation of salinity gradient at the estuary mouth (Sx|o) with the estuary 
control variables - F^ and Ft. The solid lines represent isocontours of Sx|o- 



17 




Fig. 2. Estuary Classification Diagram. The lines represent isocontours of $o- The three 
regions represent three types of estuaries: (a) light Grey- well mixed, (b) White - partially 
mixed, and (c) dark Grey- highly stratified or salt wedge. The letters denote estuaries: C - 
Columbia, J - James, Te - Tees, SW - Southampton Waterway, Ta - Tay, NM - Narrows of 
the Mersey and B - Bristol Channel. For data, see Table [1} 
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Fig. 3. Comparison between stratification at the estuary mouth obtained from theory with 
field data). 



19 



10" 



o 



10" 



10" 



■ 




1 


/^^^'o.^' O'' 13^'*^ 


1 Tamar 


a Bristol Channel 


y y^ y 9CD 


2 Puget Sound 


b Narrows of Mersey 


'^ X ^ 10 


3 Long Island Sound 


c Southampton Waterway ^ 


y yy y 


4 Narraganseft Bay 


d Tay 


yy 


5 North Ri\/er 


e Tees 


/■y^ 


6 York River 


f James 


yy 


/ o 


7 Chesapeake Bay 


g Columbia 


-X a 


/ 


8 Hudson 


h Vellar 


U -i/^ ^ 

X ^•/ 


/ 


9 Memmack 

10 Chiang Jiang 




" 


11 Ebre 




^ / 




12 Mississippi 




X v.- 




13Fraser 


'V 


- 


■ 


X 


J / 


- 


■ 


/ 6 


o / 




■ 


Xh o 


/ 


" 


y 


/^ m 


/ 


- 


X 


I ° 


s / 


Eq. (19) _ 


/ o 




/^. 


Eq. (18)with Ft = 


-i 


— Eq. (18) with Ft = 30 




o 


/ 
/ 
/ 


O Field Data from Geyer (201 0) 






• Field Data from Prandle (1 985) 


■ 




/ 


- 


; 




/ 


; 




^^B 


/ 

/ 1 


1 



10 



-4 



10" 



10" 



Fn 



Fig. 4. Comparison between our estuary classification scheme and the approximation $c 
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7F^ in (19). The grey area indicates the region where estuaries should ideally cluster. Field 



data from Geyer (2010) and Prandle (1985) are plotted for comparison with the theoretical 
predictions. 
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